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The Configuration Interaction (CI) method using a very large Laguerre orbital basis is applied 
to the calculation of the He ground state. The largest calculations included a minimum of 35 
radial orbitals for each I ranging from to 12 resulting in basis sets in excess of 400 orbitals. The 
convergence of the energy and electron-electron (^-function with respect to J (the maximum angular 
momenta of the orbitals included in the CI expansion) were investigated in detail. Extrapolations 
to the limit of infinite in angular momentum using expansions of the type AXj = Ax(J + h)~ p + 
Bx{J + i)~ P ~ + • • •! gave an energy accurate to 10 -7 Hartree and a value of (S) accurate to about 
0.5%. Improved estimates of (75) and (5), accurate to 10 -8 Hartree and 0.01% respectively, were 
obtained when extrapolations to an infinite radial basis were done prior to the determination of the 
J — > oo limit. Round-off errors were the main impediment to achieving even higher precision since 
determination of the radial and angular limits required the manipulation of very small energy and 
(8) differences. 
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I. INTRODUCTION 

Large configuration interaction (CI) calculations of the 
helium ground state are performed here in order to more 
precisely elucidate the convergence properties of the CI 
expansion for this atom. The general properties of the 
CI expansion have been known since the seminal work of 
Schwartz yj, which provided the underlying foundation 
for the later analytic and computational investigations 
HH H H H H Hi The CI expansion using single cen- 
ter orbitals is slowly convergent with respect to J, the 
maximum angular momentum of any orbital included in 
the CI expansion. In particular, the leading term to the 
energy increment is expected to behave at high J as: 

AE j = (Ey-(Ey-^jj^. (i) 

Although there have been a number of very large CI 
calculations performed on helium, all of the earlier calcu- 
lations using analytic basis sets treat the higher J con- 
tributions to the energy with less precision than the low 
J contributions. Typically, the number of radial orbitals 
for the high £ are smaller than the number of low £. The 
justification for this is that the high I partial waves make 
a smaller contribution to the energy and other expecta- 
tion values of than the low I orbitals. At first sight this 
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approach would seem reasonable for obtaining estimates 
of the total energy. However, this approach does lead 
to problems when studying the convergence properties of 
CI expansion itself. Here it is necessary to ensure that 
the successive contributions to the energy are obtained 
with the same relative accuracy and this can hardly be 
guaranteed with a radial basis that decreases in size as i 
increases. Indeed, the evidence suggests that the dimen- 
sion of the radial basis should be increased as J increases 
if the relative accuracy of the energy is to be maintained 

EJIUJ 

The convergence problems present in CI calculations 
of atomic and molecular structure are also present in a 
much severe manner in CI calculations of the positron- 
atom problem. The CI method has recently been applied 
to the study of positronic atoms (electronically stable 
states consisting of a positron bound to an atom) 0, 0, 
EHIl E E El El El E3 and also to positron-atom 
scattering states j2H |23, l22j. The attractive electron- 
positron interaction leads to the formation of a Ps cluster 
(i.e. something akin to a positronium atom) in the outer 
valence region of the atom [H, EH H3, HE} . 

The accurate representation of a Ps cluster using only 
single particle orbitals centered on the nucleus requires 
the inclusion of orbitals with much higher angular mo- 
menta than a roughly equivalent electron-only calcula- 
tion [Til IT^ . Efl W\ . In the most extreme case so far 
considered, namely e + Li, a calculation with J = 30 
was required before the energy had decreased enough 
to establish binding. Given that helium is described 
as slowly convergent [l|, one struggles to find an adjec- 
tive that could characterize the convergence properties of 



2 



positronic systems! 

The two most important expectation values for 
positronic systems are the energy, and the rate for 
electron-positron annihilation. The annihilation rate, 
which is proportional to the expectation of the electron- 
positron (5-function, has the inconvenient property that it 
is even more slowly convergent than the energy with re- 
spect to orbital angular momentum. One has successive 
increments decreasing at high J according to ITU l2S| : 



ap 7 = (iy - (iy- 1 



(J- 



(2) 



To put this in perspective, it would take a calculation 
with J « 250 to recover 99% of the PsH annihilation 
rate ^(| . In addition to the slow convergence with J, the 
(5-function operator also exhibits very slow convergence 
with respect to the radial basis [lOj . 

In the present work, large basis CI calculations of the 
He ground state are performed in order to more exactly 
understand the convergence of the CI expansion. Since 
the properties of the He ground state are known to high 
precision it is a very useful laboratory system with which 
to test methods of extrapolating the radial and partial 
wave expansions to completion. The insights obtained 
from helium should then be applicable to positronic sys- 
tems and also possibly give additional guidance about 
how to approach purely electronic systems. Besides look- 
ing at the energy, the convergence of the CI expansion of 
the electron-electron 5-function expectation value is also 
studied due to its relation with the electron-positron an- 
nihilation operator (which is also a (5-function). It should 
be noted that the (5-function operator also appears in the 
Breit-Pauli relativistic correction as the two-body Dar- 
win interaction HjJ . The present work builds on an 
earlier investigation that studied the convergence of the 
radial basis in a simplified model of the helium atom 
which only included I = orbitals |3(J. 



II. THE CI METHOD AND CONVERGENCE 
PROPERTIES 

The CI wave function in a single-center basis is a linear 
combination of anti-symmetrised two-electron states with 
the usual Clebsch-Gordan coupling coefficients, 

|*; LS) = c v M (timitjmj\LM L ) {^^\SM S ) 



1,3 



(3) 



The functions </>(r) are single electron orbitals written as 
a product of a radial function and a spherical harmonic: 



All observable quantities can be defined symbolically as 

J 



(•5) 



L=0 



where AX J is the increment to the observable that oc- 
curs when the maximum orbital angular momentum is 
increased from J — 1 to J, e.g. 



AX J = (X) J - (A)' 7 - 1 . 
Hence, one can write formally 

oo 

(A)°° = (A) J + J AX 1 

L=J+l 



(6) 



(7) 



The first term on the right hand side will be determined 
by explicit computation while the second term must be 
estimated. The problem confronting all single center cal- 
culations is that part of (A) 00 arises from terms with 
^-values that are not included in the largest explicit cal- 
culation. The two expectation values that were inves- 
tigated were that of the energy {E)°° and the electron- 
electron (5-function, {S)°° = (S(ri — r^)) 00 . For helium, 
terms with £ > 2 contribute only 0.033% of the total en- 
ergy. For purely electronic systems these higher I terms 
make a small (but slowly convergent) correction to the 
total energy and other expectation values. 

The extrapolation schemes used later in this paper 
have their basis in the work of Schwartz , Hill Q and 
Kutzelnigg and associates 0, • Analytic work indicates 
that the energy increments are given by 



AE J 



where 



-4/ 



+ 



B E 



C e 



+ 



(8) 



A E = -6tt 2 J \^{r,r,0)\ 2 r 5 dr = -0.074226 (9) 

B E = — ^ f |^(r,r,0)| 2 r 6 dr = -0.030989 (10) 

given a two-body wave function \P(ri, r 2 , (ri — r 2 )). No 
expressions for Ce have been presented. At large J, one 
expects the energy increments to be well described by 
eq. 0). 

For the 5-function one can write 



A5 J 



A s 



(r) = P(r)Y em {v) 



(4) 



where As is believed to be 

A 5 = -4vr J |f(r,r,0)|Vdr = -0.04287 (12) 

(Ottschofski and Kutzelnigg give a formula similar to 
this for the leading relativistic contribution to the en- 
ergy of two-electron atoms. We have assumed the slow 
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As/(J + |) 2 convergence is due to the two-electron Dar- 
win term). It should be noted that Gribakin and Ludlow 
|28j have also derived an expression equivalent to eq. fT^|l 
in the context of positron annihilation. The numerical 
value was taken from a variational wave function of the 
He ground state with a basis of 250 explicitly correlated 
gaussians and an energy of -2.9037243752 Hartree. 

As can be imagined there have been a number of very 
large CI calculations upon the helium ground state that 
have addressed the convergence issue 0, III IE IE IE HE] • 
These calculations can be roughly divided into two 
classes, those that represented the radial wave function 
on a grid or used piecewise polynomials 0, IE HE] > an( i 
those that describe the radial wave function as a linear 
combination of analytic basis functions 0,|EI3- The first 
systematic calculation was the seminal investigation by 
Carroll, Silverstone and Metzger (CSM) Q who used a 
piece- wise polynomial basis to construct a natural orbital 
expansion. Besides performing some very large calcula- 
tions they also estimated the completeness limit of their 
radial basis. The largest explicit calculation by CSM 
will be termed the CSM calculation while the extrapo- 
lated calculation will be denoted as CSM^. Despite their 
importance, these calculations have been largely super- 
seded by the grid-based calculation of Salomonson and 
Oster (SO) |32j and the B-spline calculation of Decleva, 
Lisini and Venuti (DLV) |5J . The SO calculation obtained 
energies, (E) J accurate to about 10 -8 Hartree by extrap- 
olating the radial basis to the variational limit. This ex- 
treme level of accuracy has not been achievable with the 
three calculations that used Slater Type Orbitals (STO) 
to represent the radial wave function [j, |E IE] • Linear de- 
pendence problems become severe as the basis set is ex- 
panded toward completeness. Indeed, recourse was made 
to very high precision (REAL*24) arithmetic in the Sims 
and Hagstrom (SH) calculation Q which is the largest 
calculation of this type so far reported. 



III. THE PRESENT CI CALCULATIONS 

The present calculations use a basis set consisting of 
Laguerre Type Orbitals (LTOs) [IE EE HE Hi- The 
LTOs of a given £ are chosen to have a common exponen- 
tial parameter which means they are automatically or- 
thogonal. Hence, the basis can be expanded toward com- 
pleteness without causing any linear dependence prob- 
lems. The CI basis can be characterized by the index 
J, the maximum orbital angular momentum of any sin- 
gle electron orbital included in the expansion of the wave 
function. It should be noted that all matrix elements 
were evaluated using gaussian quadrature even though 
the basis functions have an analytical form |15|. 

Three sets of calculations have been performed for the 
He ground state. In the first set, there were 20 LTOs 
per £ with the largest calculation including orbitals up 
to £ = 12. The LTO exponents for a given £ were the 
same and the values of the exponents were optimized in 



a quasi-perturbative fashion. The exponents for I = 0, 1 
and 2 orbitals were optimized in a CI calculation with all 
60 orbitals. The exponents for £ > 2 were optimized sep- 
arately for each £ with CI calculations that also included 
the I = 0,1,2 orbitals. Once the exponents were opti- 
mized, a sequence of calculations to give the (E) J and 
(5) J for successive J were carried out. The basis is de- 
noted the 20LTO basis and the results of the calculations 
with this basis are reported in Tabled 

The second set of calculations were much larger. Here 
there were 35 LTOs per £ with the exception of £ = and 
1 where respectively 44 and 36 LTOs were used respec- 
tively. The orbital exponents were optimized for each £ 
in a manner similar to that described above and the cal- 
culations were taken to J — 12. A total of 465 single 
electron orbitals were included in the largest calculation, 
which required the diagonalization of a hamiltonian ma- 
trix of dimension 8586. This calculation was an example 
of a very large explicit calculation. The basis is denoted 
the 35LTO basis and the results of the calculations with 
this basis are reported in Tabled 

The idea behind the third calculation was to exploit 
extrapolation techniques to estimate the variational limit 
for each partial wave. A sequence of calculations with 32, 
33, 34 and 35 LTOs per £ was done for a basis that was de- 
fined with the same exponential parameters as the 20LTO 
calculation. The number of basis functions were varied 
so that all partial waves had the same basis dimension. 
Optimizing the LTO basis for the largest radial basis has 
been shown to result in distortions in the convergence 
pattern with respect to the number of radial basis func- 
tions. This can be avoided if the basis is optimized in a 
basis that has at least 10 fewer LTOs per £ than the ac- 
tive calculation j^E] ■ The variational limit for the radial 
basis can be estimated by fitting the increments to (E) 
and (5) to the inverse series [3Cj 



AE N 
A5 N 



a-E 
as 



j\T8/2 

bs 



CE 
AT9/2 

(-S 



_/V"5/2 ^ N a/ 2 ^ N7/2 



(13) 
(14) 



It is possible to estimate the N — > oo limits for the radial 
basis once the a^, ag, b^, ■ ■ ■ coefficients have been de- 
termined. A two-term series was used for both eqs. Ijl3|) 
and i|14[l . It would have been preferable to use three- 
term series but the impact of round-off error rendered 
this impractical. The basis for this set of calculations is 
denoted the 35LTO* basis, while the basis including the 
N — > oo correction is termed the SSLTO^, basis. The en- 
ergies and expectation values for these two calculations 
are listed in Table ITT1 



IV. INVESTIGATION OF THE PARTIAL WAVE 
SEQUENCE 

The validity of these results can be tested by exami- 
nation of the energy increments of large CI calculations 
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TABLE I: Results of the present set of 20LTO and 35LTO CI calculations of He giving the energy (E) J and delta- function 
(6) J expectation values as a function of J (all energies are given in Hartree, while {8} J is in The total number of electron 
orbitals is N OT h while the LTO exponent for £ — J is listed in the A column. The results in the three {E)°° rows use inverse 
power series of different length to estimate the J — > oo extrapolation. 



20LTO 35LTO 



J 


A 


iVorb 


{E) J 






(sy 




A 






(E) J 






(sy 







4.8 


20 


-2.879 028 


507 





155 


789 


346 


8.6 


44 


-2 


879 


028 


760 





155 


766 


769 


1 


7.8 


40 


-2.900 515 


873 





128 


501 


540 


11.6 


80 


-2 


900 


516 


228 





128 


460 


082 


2 


10.1 


60 


-2.902 766 


378 





120 


923 


186 


14.4 


115 


-2 


902 


766 


823 





120 


862 


126 


3 


12.1 


80 


-2.903 320 


527 





117 


264 


315 


17.2 


150 


-2 


903 


321 


045 





117 


183 


496 


4 


14.0 


100 


-2.903 517 


973 





115 


104 


494 


19.2 


185 


-2 


903 


518 


552 





115 


004 


651 


5 


15.5 


120 


-2.903 605 


022 





113 


681 


991 


21.2 


220 


-2 


903 


605 


654 





113 


563 


078 


6 


17.1 


140 


-2.903 649 


142 





112 


676 


622 


22.8 


255 


-2 


903 


649 


820 





112 


539 


353 


7 


18.7 


160 


-2.903 673 


821 





111 


930 


245 


24.8 


290 


-2 


903 


674 


539 





111 


775 


243 


8 


20.1 


180 


-2.903 688 


677 





111 


355 


981 


26.5 


325 


-2 


903 


689 


430 





111 


183 


690 


9 


21.5 


200 


-2.903 698 


142 





110 


901 


652 


28.0 


360 


-2 


903 


698 


926 





110 


712 


681 


10 


22.9 


220 


-2.903 704 


451 





110 


534 


160 


29.5 


395 


-2 


903 


705 


263 





110 


329 


155 


11 


24.2 


240 


-2.903 708 


815 





110 


231 


642 


31.0 


430 


-2 


903 


709 


652 





110 


011 


117 


12 


25.5 


260 


-2.903 711 


927 





109 


978 


870 


32.5 


465 


-2 


903 


712 


786 





109 


743 


369 


Exact [31] 




-2.903 724 377 034 





106 


345 


371 























{E)°° and (<5)°° extrapolations 



Method 1 


-2.903 723 421 


0.106 943 


-2 


903 


724 


362 


0.106 527 


Method 2 


-2.903 723 252 


0.107 178 


-2 


903 


724 


249 


0.106 630 


Method 3 


-2.903 723 205 


0.107 334 


-2 


903 


724 


240 


0.106 698 



35LTO^ basis sets, along with the SH, CSM^, SOoo and 
DLV calculations. These same sets of data are also pre- 
sented as energy differences between consecutive calcu- 
lations AE J . The energies of the SH calculation, which 
used the even-tempered STO basis, are consistently the 
worst, and are 3 x 10 -7 Hartree larger than the 35LTO 
calculation at J = 12. Even though CSMqo does attempt 
to achieve the variational limit for each J, in reality it is 
only about as good as the 35LTO calculation. Indeed for 
J > 7, the CSMqo values of AE J were smaller than those 
of the 35LTO calculation. The present 35LTO^, SOoo 
and DLV calculations are in agreement to 10~ 8 Hartree 
(or better) for J < 4. This is expected, since all three 
calculations are large and use extrapolation techniques 
to achieve the variational limit. The energy difference 
between the 35LTO and 35LTOoo energies gives an indi- 
cation of the incompleteness of the 35LTO basis and by 
,/ = 12 the difference is 0.96 x 1CP 7 Hartree. 

The good agreement between the 35LTO^, SOoo and 
DLV energies is not present for J > 5. Although the 
35LTO^, and SOoo energies generally agree at the level 
of 1CP 8 Hartree, it is seen the DLV AE J increments are 
larger than these two other calculations. For example, 
DLV give AE 5 = 8.737 x 1(T 5 Hartree which is about 
2 x 1(T 7 Hartree larger than the 35LTO^ and SOoo in- 
crements. It has also been noticed that DLV do overstate 
the accuracy of their calculation, they assert an accuracy 
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FIG. 1: The exponents pe as a function of J for the different 
CI calculations of the He ground state energy as listed in 
Table ITTT1 



of helium. Besides the present calculations, data from a 
number of previous CI calculations have been used. 

Table ITO! gives the energies of the present 35LTO and 



5 



TABLE II: Results of the 35LTO* and 35LTO£o CI calculations of He for the (E) 3 and (8) J expectation values as a function of 
J (all energies are given in Hartree, while (<5) J is in a§). The total number of electron orbitals is A^b while the LTO exponent 
for £ = J is listed in the A column. The results in the three (i?) 00 rows use inverse power series of different lengths to estimate 
the J — * oo extrapolation. 



35LTO* 35LTO* 



J 


A 


iVorb 


(E) J 










(E) J 






(S) J 





4.8 


35 


-2.879 


028 716 


0. 


155 


774 


273 


-2. 


.879 


028 


766 


0. 


155 


763 804 


1 


7.8 


70 


-2.900 


516 172 


0. 


128 


472 


171 


-2. 


900 


516 


246 


0. 


128 


451 020 


2 


10.1 


105 


-2.902 


766 757 


0. 


120 


878 


722 


-2. 


.902 


766 


852 


0. 


120 


845 876 


3 


12.1 


140 


-2.903 


320 971 


0. 


117 


204 


759 


-2. 


903 


321 


084 


0. 


117 


159 843 


4 


14.0 


175 


-2.903 


518 472 


0. 


115 


030 


202 


-2. 


903 


518 


601 


0. 


114 


973 165 


5 


15.5 


210 


-2.903 


605 568 


0. 


113 


593 


010 


-2. 


903 


605 


710 


0. 


113 


523 543 


6 


17.1 


245 


-2.903 


649 729 


0. 


112 


573 


360 


-2. 


903 


649 


884 


0. 


112 


491 488 


7 


18.7 


280 


-2.903 


674 443 


0. 


111 


813 


216 


-2. 


903 


674 


609 


0. 


111 


719 074 


8 


20.1 


315 


-2.903 


689 330 


0. 


111 


225 


558 


-2. 


903 


689 


506 


0. 


111 


119 165 


9 


21.5 


350 


-2.903 


698 823 


0. 


110 


758 


273 


-2. 


903 


699 


007 


0. 


110 


639 719 


10 


22.9 


385 


-2.903 


705 157 


0. 


110 


378 


300 


-2. 


903 


705 


349 


0. 


110 


247 727 


11 


24.2 


420 


-2.903 


709 543 


0. 


110 


063 


722 


-2. 


903 


709 


741 


0. 


109 


921 236 


12 


25.5 


455 


-2.903 


712 675 


0. 


109 


799 


333 


-2. 


903 


712 


882 


0. 


109 


645 079 










(E)°° and (S)°° 


Extrapolations 
















Method 1 






-2.903 


724 243 




0.106 881 


-2. 


903 


724 


476 




0.106 328 


Method 2 






-2.903 


724 123 




0.106 757 


-2. 


903 


724 


378 




0.106 341 


Method 3 






-2.903 


724 109 




0.106 847 


-2. 


.903 


724 


384 




0.106 334 



of 7.8 x 10 8 Hartree. However, this accuracy is based on 
a calculation which gives {E) 3 = -2.903319811 Hartree 
(2nd column of Table IV of 0), and this energy is in 
error by 1.3 x 10~ 6 Hartree! 



A. Scrutiny of the partial wave increments 

A useful way to scrutinize the partial wave series is to 
assume a power law decay of the form 



AX 



(15) 



and determine the value of p for a succession of three 
(X) J values using 



, = ,„(» U' J 



\ AX J 



J 



(16) 



The exponent derived from the energy increments is pe 
while the exponent derived from the ^-function incre- 
ments is p/>. One expects pe — > 4 Jl| and ps — > 2 as 
J — > oo m agreement with eqs. JSJ and Qllfl. 

The values of pe for the He energies presented in Ta- 
ble ^Q] are plotted in Figure Q] as a function of J. One of 
the noticeable features of Figure ^ are the irregularities 
in some of the calculations, e.g. the SO^, CSMqo and 
351710^ calculations. The fluctuations in the present 



351710^, curve are due to the impact of round-off error on 
the radial extrapolations. The determination of the coef- 
ficients in eq. (|13fl involves the subtraction of the energies 
for calculations that differ by a single LTO. The resulting 
energy differences are very small and therefore are sus- 
ceptible to the essentially random errors resulting from 
round-off that gradually accumulate during the course of 
the computations. The irregularities in the CSMqo and 
SOqo curves are a consequence of the number of digits at 
which the energies were published @,H2- Plots of pe vs 
J were examined (but not plotted in Figure for some 
calculations 0, Q that used an STO basis set. These 
plots of pe showed much larger fluctuations than any of 
the calculations depicted in Figure ^ 

The smaller 20LTO and CSMqo calculations had plots 
of pe vs J that tended to level out at pe ~ 4.05. Indeed, 
the tendency for 20LTO trajectory to curve up indicates 
that the successive AE J increments are decreasing too 
quickly at the higher J values. The larger SOoo, DLV, 
35LTO and 35LT0£ o calculations have pe versus J tra- 
jectories that steadily decrease with increasing J and ap- 
pear to be approaching the expected limit of pe = 4 
although this is obscured somewhat for the SOoo and 
35LTOoo curves. It will be demonstrated later that the 
behavior of the 20LTO and CSMqq curves is due to slower 
convergence of the radial basis at high I. 

The tendency for pe to approach the limiting value of 
4 from above is a consequence of the fact that the Ae 
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TABLE III: Comparison of different CI calculations of the He atom ground state energy (E) J as a function of J. To 
aid interpretation, the AE J energy differences are also tabulated. The energies of the SOoo, CSMoo and 35LTOJo are the 
estimated energies in an infinite radial basis. The 35LTO^ (E) J and AE J are the smoothed values. The DLV energies for 
J > 5 are obtained by adding the AE J from Table IV of Ref. p| to their estimate of (Ef. 



J 35LTO 35LTO^, SH [8] CSM^ [2] SOoo [32] DLV [5] 



(E)° 


-2 


.879 


028 


760 


-2. 


.879 


028 


766 


-2 


.879 


028 


757 


-2.879 028 765 


-2 


.879 


028 


77 


-2.879 028 767 


(Ef 


-2 


.900 


516 


228 


-2 


.900 


516 


246 


-2 


.900 


516 


220 


-2.900 516 220 


-2 


.900 


516 


25 


-2.900 516 245 


(Ef 


-2 


.902 


766 


823 


-2 


.902 


766 


853 


-2 


.902 


766 


805 


-2.902 766 822 


-2 
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and Be coefficients of eq. JHJ have the same sign. The 
coefficients Ae and Be are derived from second and third 
order perturbation theory respectively 0,0 and have the 
same sign due to repulsive nature of the electron-electron 
interaction. One surmises that a mixed electron-positron 
system, with its attractive electron-positron interaction 
should have be — > 4 from below, and this is indeed the 

case 13 US E3, m. 

The incremental exponent for the ^-function, p$, is 
shown in Figure for the 20LTO, 35LTO and 35LTO^ 
basis sets. It should be noted that the values of ps were 
sensitive to the precision of the calculation. Originally, 
the diagonalization of the Hamiltonian was performed us- 
ing the Davidson algorithm |35| . However, this method 
could not give (5-function expectation values to better 
than 8 significant figures (irrespective of the convergence 
tolerance for the energy) . This lead to noticeable fluctu- 
ations in the p$ versus J plot. The diagonalization was 
subsequently performed using the EISPACK libraries, re- 



ducing the size of the fluctuations. The trajectories of the 
20LTO and 35LTO calculations do not appear to be ap- 
proaching the ps — > 2 limit as J — > oo. The 20LTO curve 
has a ps trajectory that diverges from 2 for J > 6 while 
the 35LTO curve diverges from 2 for J > 7. The 20LTO 
basis gives ps = 2.155 at J — 12 while the 35LTO basis 
gives ps — 2.064 at this J value. 

However, the plot of ps based on the 351^0^ se- 
quence does exhibit the correct qualitative behavior as 
J increases. The N — > oo corrections have a larger im- 
pact on (5) J than on (E) J since the former converges as 
0(N- 5 / 2 ) while the latter converges as 0(N- 7 / 2 ). The 
J = 12 value of ps was 2.008. 

The behavior exhibited in Figures ^ an d El can be at- 
tributed to the convergence of the radial basis. A larger 
radial basis is required to predict successive AE J incre- 
ments as J increases. The ratios 
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FIG. 2: The exponents pg as a function of J for the LTO 
calculations of the He ground state (<5) . 
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FIG. 3: The 35LTO:20LTO and 35LTO^ :20LTO ratios of the 
increments to {E) J and {8) J (refer to eq. 1171 1 as a function 
of J for the calculations of the He ground state. 



gives a measure of the relative impact of the 20LTO, 
35LTO and 351^0^ basis sets to a J energy increment. 
One can define a similar ratio, i?/, for the ^-function A5 J 
increments. 

The behavior of these ratios versus J are depicted in 
FigureOl Both R J E and Rj increase steadily with J. The 
slower convergence of the energy at higher J is also ev- 
ident in Table I of and explicit comment about this 
point has been made previously [jj. The (S) J expec- 
tation value is much more sensitive to the increase in 
the dimension of the radial basis and there was a 9.3% 



FIG. 4: The extrapolated J — > oo limit for the He ground 
state energy {E} 00 using three different methods to complete 
the partial wave series. The input (-E) J upon which the ex- 
trapolations were based were those of the 35LTO and the 
smoothed 35LTO^, calculations. The horizontal line shows 
the exact helium energy [3l|. 



increase in AS 12 between the 20LTO and 35LTO^ cal- 
culations. The corresponding increase in the AE 12 was 
only 0.8%. This extra sensitivity of (S) J is something we 
have noticed in calculations of positron annihilation rates 
in positron-atom systems [la . 1171 l2ll l22j even though ex- 
plicit mention of this point has not been made. 



B. Smoothing of the 35LTOJO energies 

It is apparent from Figures ^ an d El that including 
the radial extrapolations has resulted in irregularities ap- 
pearing in the SdUTO^ energy sequence. These irregu- 
larities are of order 10~ 9 Hartree at J = 12 and should be 
removed before the J — * oo corrections are determined. 

Examination of Figure |3] suggested that (R E — 1) oc 

(AE J \ 
. "T TOSc - 1 to 
20LTO J 

a G + H x J s functional form was performed over the 
J s [5,12] interval. An adjustment to the (E) J energy 
sequence was made once G, H and s were fixed. The 
values of s was approximately s w 2.85. 

The 35LTO^ values of (E) J and AE J given in Ta- 
ble [ffi] are those of the smoothed energy sequence. The 
largest change to any of the energies was 2 x 10~ 9 Hartree. 



C. Extrapolation of the partial wave series 

One of the major aims of this paper was to determine 
whether it is possible to extract the J — > oo limit from a 
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finite sequence of calculations. To this end, fits of inverse 
power series of different lengths are made to sequences 
of (X) J data, and then those inverse power series are 
summed to infinity. 

Equations (JSJ and Ullfl are the working equations. 
Fits are performed retaining just the leading order term 
(Method 1), the first two terms (Method 2), and the first 
three terms (Method 3) of these series. The fits of these 
equations use the minimum information necessary. So, 
method 1, which only retains the first Ae term of eq. JSJl, 
only requires two successive values of {E) J to determine 
Ae- Three successive values of (X) J are used to deter- 
mine Ax and Bx when the two leading terms of eq. (jHJ 
or eq. (fTTf) are used. Four successive values of (X) J are 
used to determine Ax , Bx and Cx when the three lead- 
ing terms of eq. (JHJ or eq. are used. The fits to 
determine Ax and/or Bx and/or Cx can be done to dif- 
ferent sequences of J values as a self-consistency to check 
that the two-term fits to the J = 8, 9, 10 or J = 10, 11, 12 
sets of (X) J give answers that are numerically close. 

Once the coefficients of the inverse power series have 
been determined, the J —> oo contribution is determined 
by a two-step procedure. Firstly, the series JHJ and 111|) 
are summed explicitly up to J+200. The remainder from 
> J + 201 is determined using the approximate result: 



.=.7+1 



1 



+ (p-l)(J + l)P- 



(18) 



Eq. I|18|) can be regarded as an approximation to the 
+ h)~ 2 integral using the mid-point rule. This 
approximation is accurate to 0.1% for p = 2 and J = 7. 

Figures 0] and [S] show the behavior of the extrapolated 
E and 5-function for the three different extrapolations 
as a function of J. Tables [i] and |n] gives estimates of 
(E)°° and (5)°° using the calculated values at the largest 
possible J values to determine the J — > oo corrections. 

Figure 01 shows that the quality of the 35LTO energy 
extrapolation using method 1 is inferior to methods 2 
and 3 which give (E)°° energies in agreement which each 
other at the 10~ 9 Hartree level for J > 8. However, 
using the 35LTO energies in conjunction with methods 
2 and 3 gives {E)°° values that are too large by about 
10~ 7 Hartree. This is a consequence of using a large but 
not quite complete radial basis. The use of the 351^0^ 
energies results in an energy limit that is an order of mag- 
nitude more precise than those of the 35LTO basis. Using 
method 2 for the J = 10, 11, 12 35LTO^ energies gave 
(E)°° = -2.903 724 378 Hartree, an energy that is in error 
by 10 -9 Hartree. The J — > oo corrections were only made 
using method 2 since the more sophisticated method 3 is 
more sensitive to the imperfections of the smoothed data 
sets. The smoothed energy sequence is probably not a 
perfect reproduction of the actual sequence and there is 
a tendency for the (E) 00 limit to be more negative than 
the exact energy. The method 3 estimate of (E)°°, at 
J = 12, namely -2.903 724 384 Hartree, is about 10" 8 
more negative than the exact energy. A similar level 
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FIG. 5: The extrapolated J — > oo limit for the He ground 
state (5)°° using three different methods to complete the par- 
tial wave series. The horizontal line shows the value of Drake 

m. 



of accuracy was achieved in the earlier SOoo calculation, 
their estimate of the energy in the J tooo limit was -2.903 
724 39 Hartree jH- 

The difficulties in obtaining sub 0.1% accuracy in (5)°° 
for the 35LTO sequence are readily apparent from Fig- 
ure As one increases J, the estimates of (<5)°° also 
increase and the discrepancy with the accurate value of 
Drake [3lJ gets larger. The ultimate accuracy achievable 
for the 35LTO basis is between 0.1 and 0.5%. The ap- 
parent superiority of method 1 for the 35LTO basis arises 
because errors resulting from a finite dimension radial ba- 
sis act to partially cancel errors that arise from this least 
sophisticated J — > oo extrapolation. 

However, usage of the SSLTO^ sequence permitted a 
much more accurate extrapolation to the J — > oo limit. 
The 35LTO^ basis gives estimates of (S)°° that are two 
orders of magnitude more precise. The method 2 ex- 
trapolation was only 6 x 10~ 6 ajj larger than the exact 
value [U • While the radial extrapolations did introduce 
fluctuations into the (5) J values, the relative size of the 
individual AS J increments were much larger than the 
AE J increments and thus did not lead to fluctuations in 
(S) 00 as long as method 2 was used for the angular ex- 
trapolations. However, the use of method 3 did result in 



fluctuations of order 
in Figure [SJ 



10- 



in (5) 00 so is not depicted 



D. The coefficients of the inverse power series 

The coefficients of the asymptotic forms, © and l|ll|) 
are known a-priori from eqs. 10, (|10|) and (|12f> . Esti- 
mates of these parameters are also obtained during the 
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FIG. 6: The value of Ae as extracted from sequences of 
(E) J data. The horizontal line shows the value of from eq. ©, 
namely -0.074226. Estimates of Ae from Method 2 are drawn 
with dashed lines while estimates of Ae from Method 3 are 
drawn with solid lines. The DLV data analyzed here was 
solely from their Table IV of Ref.0, since this avoided the 
discontinuity at J = 4 and gave a smooth curve. 



fit of the inverse power series to a set of (E) J or (S) J . An 
ideal consistency check would be estimates of Ae and As 
that steadily approached -0.074226 and -0.04287 as J in- 
creased and as the number of terms included in eqs. (JSJ 
and Hll(l increased. Unfortunately, this has not yet been 
achieved. The least squares analysis of the CSMqo ener- 
gies gave A E = -0.0740 and B E = 0.031 0. However, 
this value of Ae is only achieved when using (E) J for 
J £ [5,8]. The very large calculations of DLV reported 
A E = 0.07415 and B E = 0.0317 @. However, a cursory 
examination of Figure which depicts values of Ae ob- 
tained from three successive (E) J energies demonstrates 
that their value of Ae is not converging to -0.074226 
with increasing J. Applying the more sophisticated 3- 
term inverse power series to the DLV energies leads to 
an Ae that exhibits a 4% variation between J = 6 and 
J = 13. 

Figure also shows the variation in Ae when fitting 
the 35LTO and 35LTO|^ energy sequences to eq. (0. 
Fits were performed with both methods 2 and 3 for the 
35LTO energy sets, and only with method 2 (for rea- 
sons discussed earlier) to the 351^0^ sequence. The 
Ae coefficients for a given method are computed using 
the minimum range of J values that permitted the unique 
determination of the coefficients. 

Application of method 2 to the 35LTO data reveals 
that Ae achieves a minimum value of Ae = -0.07413 
at J = 7 before increasing at larger J. Application of 
method 3 results in values of Ae that are clearly not 
approaching the correct value. This should be expected 



FIG. 7: The value of As as extracted from sequences of (5} J 
data. The horizontal line shows the value of eq. 1121 . namely 
A s = -0.04287. 



since it has been demonstrated that the 35LTO AE J 
are increasingly underestimated as J increases. It would 
therefore be hoped that values of A J E extracted from 
the 35LTO^, would show better convergence to the ex- 
pected limit as J increases. This expectation has only 
been partly realized, there are indications that Ae may 
be converging to the correct value, but the application 
of smoothing has probably introduced a systematic bias 
that resulted in a tendency to overestimate the magni- 
tude of A E - 

Figure [7| shows the values As as obtained from the 
35LTO basis using Methods 1, 2 and 3 as a function of J. 
None of the calculations using the 35LTO basis resulted 
in an As vs J curve that approached the correct value 
as J increased. This is another manifestation of the very 
slow convergence of (5) J with respect to the dimension of 
the radial basis set. There was a significant improvement 
when As was extracted from the 35LT0£ o sequence using 
method 2. In this case, As does appear to be converging 
to the expected value of -0.04287 and at J = 12 one 
obtains A s = -0.04282 

The small irregularities in the i^LTO^ (5) J sequence 
resulted in irregularities in As when using the more so- 
phisticated and sensitive method 3 fit; so this was not 
depicted in Figure [7| It is also worth noting that As was 
also subject to irregularities of ±2% when the David- 
son method was used to diagonalize the hamiltonian and 
generate the ground state wave function. 



V. SUMMARY AND CONCLUSIONS 

Results of a set of very large CI calculations of the He 
ground state have been presented. The largest explicit CI 
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calculation reported here with a minimum of 35 LTO's 
per i gave an energy that was accurate to 1.2 x 10~ 5 
Hartree. Including energy corrections obtained from the 
2- and 3-term inverse power series in J resulted in a He 
ground state energy that was accurate to « 1 x 10~ 7 
Hartree. Improved accuracy required the use of extrap- 
olations in the radial basis set to get an estimate of the 
variational limit for (E) J . This permitted the He ground 
state energy to be predicted to better than 10~ 8 Hartree, 
an improvement of a factor of 1000 over the largest ex- 
plicit calculation. The main impediments to more refined 
predictions of the He ground state are those due to round- 
off errors. Estimating the coefficients of the inverse power 
series involves manipulating very small energy differences 
which will be sensitive to round-off errors. The fluctu- 
ations in the radial extrapolation were about 2 x 10~ 9 
Hartree at J = 12. While this in itself is not that bad, 
these fluctuations are magnified by an order of magni- 
tude when the angular momentum extrapolation is then 
done. The impact of the fluctuations was somewhat mit- 
igated by the introduction of a smoothing procedure, at 
the cost of introducing a small systematic error. 

The prediction of the electron-electron (5-function was 
considerably more difficult due to the 0(L + i) 2 conver- 
gence. In this case, the explicit calculation was accurate 
to 3% at J — 12. Application of the inverse power series 
(Method 2) to include higher J contributions improved 
the accuracy to 0.3%. The main reason for the low accu- 
racy was the slow convergence with respect to the num- 
ber of radial basis functions. The relative accuracy of 
successive AS J increments decreases as J increases if the 
number of radial basis functions per I is kept the same. 
Once again, extrapolating the radial basis to the varia- 
tional limit lead to an improved prediction of (S). The 
best CI estimate of (5) = 0.106341 was within 0.01% 
of the close to exact variational estimate |3lJ . The ex- 
trapolations of ((5) were less susceptible to round-off error 
simply because the A5 J increments were larger. 



While the use of extrapolations did improve the qual- 
ity of the calculation, the full potential of the method 
has not been realized due to round-off error. The radial 
matrix elements are evaluated with gaussian quadratures 
and the achievable precision for the larger calculations is 
about 10 -12 Hartree. This accuracy could be improved 
by the either the development of a convenient analytic 
form for the electron-electron matrix elements or the us- 
age of quadruple precision arithmetic. This would then 
permit the use of inverse power series with more terms 
leading to improved radial and angular extrapolations. 
For example, an accuracy of 10~ 14 Hartree was achiev- 
able for a CI calculation restricted to I — orbitals [3(| • 

These results have implications for the prediction of 
the annihilation rate of positronic atoms from single- 
center CI type calculations ^(|- Some sort of extrap- 
olation in J is needed to determine the energy and more 
particularly the annihilation rate. One way to minimize 
the impact of the extrapolation in J is to run the calcu- 
lation to the highest possible angular momentum. How- 
ever, the high J parts of the annihilation rate will tend 
to be increasingly underestimated as J increases unless 
accurate estimates of the radial variational limit can be 
made. Since this can now be achieved for a Laguerre 
basis [33 , it eminently conceivable that estimates of the 
annihilation rate at better than 0.1% accuracies will be 
achievable for single-center basis sets. 
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